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Abstract. This note attempts to interpret some of the recent findings about a downtrend 
in the mean azimuthal velocity of low [ff/Fe] thin disc sta rs with increasing metallicity . 
The presence of such a trend was predicted in the model of lSchonrich & Binnevl (2009a), 
albeit with a slightly steeper slope. We show that in a simple picture a Galactic disc with- 
out mixing in angular momenta would display an exceedingly steep trend, while in the 
case of complete mixing of all stars the trend has to vanish. The difference between model 
and observational data can hence be interpreted as the consequence of the radial abun- 
dance gradient in the model being too high resulting in an underestimate of the migration 
strength. We shortly discuss the value of eccentricity distributions in constraining struc- 
ture and history of the Galactic disc. 

1 About the trend in azimuthal velocities 

Among the various links between chemistry and kinematics in the solar neighbourhood is the trend in 
mean azimuthal velocity of stars with increasing metallicity for objects belonging che mically to the 
thin disc. The observational evidence for this trend can be traced back to lHavwoodl (I2008I) . who found 
that metal-poor st ars at low [g/Fe] in local sam ples have large angular momentum. 

As laid out in Schonrich & Binnev (2009b) the low [ g/Fel stars cons t itute th e younger part of the 
Galactic disc (cf. Matte ucci & Brocatolll990l) . The lSchonrich & Binnevl d2009al) model cannot make 
a firm prediction about such trends for the high [a/Fe] stars that formed at early epochs. It is still too 
uncertain what the radial abundance gradient and the radial behaviour of star formation approximately 
looked like. In particular either an inverse radial metallicity gradient at early times or different setups 
of inside-out formation can give rise to an inverse (i.e. positive) gradient of mean azimuthal velocity 
with [Fe/H] for the older populations. We reserve the discussion of such possibilities to a paper in 
progress. Yet for the younger, low [q/F eJ stars very firm pred ictions are available (for a discussion 
of these effects in an N-Body model see Loeb man et all |20TH) . Observations imply a very flat age- 
metallicity relation during the past couple of Gyrs. Hence the links between chemistry and kinematics 
are dominated by the radial abundance gradient in the interstellar medium: With increasing age stellar 
populations acquire random energy, which brings them on more eccentric orbits letting them oscillate 
further away from the circular orbits defined by their angular momentum. We call this "blurring". As 
angular momentum is conserved in an axisymmetric potential, we see in the solar neighbourhood stars 
from inner radii at low azimuthal velocities and vice versa, so that the radial abundance gradient gives 
rise to a significant downtrend of mean azimuthal velocity with metallicity. 

Radial migration shuffles stars in angular momentum, a process called "churning", but the popula- 
tions should still contain some residual information on where they were born: as stars from the inner 
Galaxy get scattered to larger angular momenta from where they can visit the Solar neighbourhood, 
the metal-rich population is still preferentially concentrated towards the inner Galactocentric radii, 
while the metal-poor outer disc populations are preferentially at large angular momentum. The expla- 
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nation in Scho nrich & Binnevl d20 09b) focused on the fact that radial migration is largely responsible 
for us being able to find those migrated objects in large numbers at all and on that we can observe 
this gradient in contrast to classical expectation where the increasing asymmetric drift with age should 
have given an inverse gradient from an sloping age-metallicity relation. Unfortunately it evoked the 
widespread misconception that radial migration is responsible for the downtrend of azimuthal veloc- 
ities with metallicity. This is not true, because in fact radial migration tends to weaken the observed 
trend. 

We can illustrate this by t wo extreme cases: Assume a radial abundance gradient of -0.06 dex/ kpc 
(see iLuck & Lambert], 1201 ll) in the young disc. In the first case we assume that stars do not change 
their initial angular momentum and that there is insignificant scatter of stellar abundances at birth at a 
certain Galactocentric radius R g . We can hence write for the young population: 

[Fe/H](^) = [Fe/H] + ^^(R g - Rq) (D 

UKg 

where Rq = 8.2 and [Fe/H] denote the Solar galactocentric radius and the local metallicity. Assuming 
a disc with constant circular velocity of V c = 233 kms~' we can then replace the galactocentric radius 
R g using angular momentum conservation and resolve the linear relationship for the local azimuthal 
velocity V: 

R q l Ra-Ro\ I [Fe/H] - [Fe/H] dR q \ 
V ([Fe/H]) = V c -1 = V c 1 + - = V e 1 + 1 ' 1 - 1 ° * (2) 



Ro \ Rq I \ Rq d[Fs/U] 

and for the chosen parameters we obtain the slope in (on the quite meagre baseline in metallicity 
that would be observable): 

dV VcdRn , 

^^ = " 474kfflS /d6X - (3) 

The res ult would be about dV,p / c/[Fe/H] ~ -300 km s _1 / dex for the values used in lSchonrich & Binnevl 
d2009al) . largely due to the steeper metallicity gradient assumed there. Measurement errors in metal- 
licity would blur out the baseline and reduce the measured gradient, but this will not remove the stark 
contrast to what is observed in the data, as well as contributions from age-dependent evolution of the 
Galactic disc can hardly remove all difference on the super-solar metallicities. Radial migration with 
its redistribution in angular momentum, however, does a huge change to this relation: If we assume 
the other extreme case of complete radial mixing, the gradient will vanish, as the distribution of pop- 
ulations does not depend any more on where the stars have been born. So we see that in general radial 
migration reduces the gradient and is both detectable by the expansion of the baseline in metallicity 
and by the reduced gradient. The presence of a weak velocity-metallicity correlation in data confirms 
that there is considerable radial migration, but the mixing is not com plete. 

Let us examine those trends on a real data set. From the model of lSchonrich & Binnevl (l2009al) we 
expect that if we do not discard alpha rich stars the thick disc velocities will moderate any trend in 
the mean azimuthal velocities with mean metallicities. In the model they balance the young disc trend 
at intermediate metallicities, because the downtrend for the younger low [ff/Fe] stars is compensated 
for by the increasing number of high asymmetric drift, high [ff/Fe] stars of the chemical thick disc 
towards lower metallicities. In the lower metallicity regime ([Fe/H] < -1.0) we even encounter a 
mildly positive slope in the model, because of some peculiarities in early disc formation that will be 
discussed in an upcoming paper. In this light it is essential to get a clean cut in the metallicity plane, 
removing the high [ff/Fe] population. The following exercise was first presented at the IAU assembly, 
but concerning the need for re-interpretation it seems appropriate to show it again. 

Fig. [D shows the m etallicity plane as defined by [Mg/Fe] as alpha element against [Fe/H] in the 
Borkova & Marsakov (2005) sample. This sample is a homo genised compilation of high-resolution 
spectroscopic data from several studies of local stars (e.g. [Bensbv et al., 2005; Fuhrmanrl 120041 : 
Red dv et alll2003l l2006h . One should keep in mind a major caveat against over-interpretation of any 
data we use here: Almost all local high-resolution data have strong kinematic selection biases for in- 
creasing or decreasing the number of selected thin/thick disc stars. It should also be kept in mind that 
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Fig. 1. The metallicity plane in [Mg/Fe] versus [Fe/H] of the Borkova et al. sample. The depopulated region in 
the data advises a cut somewhere between [Mg/Fe] = 0.14 and [Mg/Fe] =0.18 where we set the cut. 



because of the strong correlations between single velocity components and also between metallicities 
and kinematics, any such sample is uncontrollably biased both in metallicities and in velocity space. 
Yet, it is still interesting to look at these trends, which are confirmed on GCS data in these proceed- 
ings by L. Casagrande. As advocated by the depopulated region at intermediate [Mg/Fe] we cut the 
sample to keep only objects with [Mg/Fe] < 0.18, to get a chemical thin disc selection. In light of 
the observational errors there might be some residual contamination, which seems to be confirmed 
when we measure the gradient of metallicity against angular momentum: Since the local azimuthal 
velocities can be translated directly into the angular momentum of those stars and let us this way es- 
timate the stellar guiding centre radii, a local sample can be used as a good indicator of the radial 
abundance gradient. For our cut we estimate d[Fe/H]/dR = -0.04dex/kpc a value that rises up to 
-0. 1 dex/ kpc depending on how much of the low angular momentum regime we remove from the 
sample and how far down in [Mg/Fe] we move the cut (the radial gradient increases mildly for stricter 
selections which indicates that we are tossing out s ome residual old s tars). Interestingly the gradient 
for the cleaner selections is a bit higher than in the iLuck & Lambertl d201 lb data, but this should not 
be taken too seriously in the light of the diverse biases in our sample. 

The trend of mean heliocentric V velocity against metallicity in the low [Mg/Fe] subsample is 
examined in Fig. [2] A linear fit (green line) yields a highly significant slope of dV,p/d[FefH] = 
(-31.5+ 4.2) km s" 1 / dex a nd is plotted onto the data. The observed trend is a bit higher than what was 



found bylLee et al. (1201 II) - this can most likely be traced to the larger [a/Fe] errors of the SEGUE 
(lYannv et all 200S ) sample with its low resolution spectra that likely results in larger contamination 
of the low [a/Fe ] part w ith high [a/Fe] stars. Surprisingly our finding is not consistent with the result 
of lNavarro et alJ d201lh on essentially the same data. The same result as with the linear fit is seen from 
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Fig. 2. Kinematics versus metallicities in the Borkova dataset in the subsample of stars with [Mg/Fe] < 0. 18 (red 
data points). The blue line gives the linear fit of the sample compared to the binned heliocentric mean velocities 
shown in green. The purple line shows the velocity dispersions. 



the binned means, where we divided the sample into 0. 1 dex wide subsets. At the same time we plot the 
dispersions for each bin (purple line). We can compare this with the qualitativ e predictions of the dif- 
ferent chemical evolution models: In classical chemical evolution models (e.g. lChiappini et aTll200ll) 
the thin disc density ridge is created by the local population running along this ridge from low to high 
metallicities. This implies a significant age trend with metallicity and hence via the age-dispersion 
relation the most metal poor stars should show the highest velocity dispersions with a clear downtrend 
towards the metal rich objects. This is not observed, while the data are consistent with the behaviour 
in the radial migration models that favour dispersions being rather flat with metallicity at low [ar/Fe]. 
Similarly the mean V velocity trend is hard to explain in the classical framework, but straight forward 
in the framework of radial migration models. Yet there is a little problem: The trend of V velocities 
with metallicity is considerably higher in the model (predicting up to 50kms -1 / dex) than it is in the 
data. 

When we look at th e above implication s of the velocity gradient on metallicities, it becomes clear 
what went wrong in ISchonrich & Binnevl (I2009al) . When we set the radial abundance gradient in 
the model, we had to rely on older abundance gradient data that now seem to hav e indicated a too 
steep abundance gradient. At least it is significantly steeper than the value derived bv lLuck & Lambertl 
(1201 ll) . The overestimated gradient places one of the main findings of the first paper on firm grounds, 
namely that radial migration is necessary for explaining the local metallicity distribution. Yet, it ap- 
pears that with the too steep gradient we underestimated the need for stellar radial migration: The 
lower the gradient, the farther stars need to migrate in order to yield the same width of the local 
metallicity distribution. Strengthened migration should then result in a shallower slope of velocities 
with metallicity. As a side effect the stronger migration would result in a locally stronger and larger 
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Fig. 3. The eccentricity distribution from Schonrich & Binnev ( 2009a) at different altitudes above the plane at 
Solar Galactocentric radius. We used simple population masses, i.e. no selection function was applied. 



scale-height thick disc. The standard model from I Schonrich & Binnevl <|2009a) was already rather on 
the upper edge on how strong the thick disc may be. Again thi ngs fit together as on the dec ision be- 
tween vertical energy and vertical action as conserved quantity Schonrich & Binnev (2009a) decided 
for vertical energy. Solway et al. (in prep) and also the recent paper of iBird et al.l d201 lb showed that 
vertical action is rather conserved than vertical energy. This will deliver some reduction in thick disc 
scale heights that can balance the stronger thick disc arising from increased radial migration. 



An impediment to modelling has so far been that despite its advantages th e classic adiabatic ap- 
proximation and also adiabatic modelling of stellar populations as put forward in lBinnevI d2010h violate 
total energy conservation. The main effect of this short-coming is a significant underestimate of the 
Galactocentric radii populated by an orbit or in other words an under-prediction of inner disc stars in 
the intermediate and out er disc regions. Two solut ions to this problem have been pr oposed: angular mo- 
mentu m correction bv lBinnev & McMillan! d201 ll) and the adiabatic potential of lSchonrich & Binnevl 
(1201 ll) . which directly restores total energy conservation by correcting the horizontal potential for the 
energy removed from the vertical motion. This solution pushes the high vertical energy orbits back to 
the outer disc, compensating locally for some of the differen ce between isotherma l appro ximations and 
the adiabatic approximation. We are currently updating the Schonr ich & Binnevl (l2009al) model to in- 
clude the adiabatic potential. A definitive treatment of the ([Fe/H], V$) relation must await completion 
of this project. 
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2 About eccentricity distributions 

In recent studies starting with LSales et all (l2009h and e.g. later iLee et al ] d201lh the comparison of 
eccentricity distributions has become a frequently used tool to assess pos sible scenarios for the history 
of the Galactic disc. Fig. [3] shows the eccentricity distribution from the Sch onrich &Binnev (2009a) 
model without change of any parameter and apply ing no selection f unction, i.e. using the population 
masses as weight. The result looks quite similar to lLee et ail d201 ll) . yet we do not want to book this 
as a win, as we would like to advocate against the use of eccentricity distributions for several reasons: 
First it is counter-intuitive to fold the wealth of information we have into a single quantity that hides 
essential information like the angular momentum distribution of the disc. As another point one needs to 
perform a real propagation of distance, proper motion and line-of-sight velocity errors that to date does 
not appear to have been fully carried out. This would be computationally expensive, because in contrast 
to simple velocity space errors have to be calculated on the orbital model from which the eccentricities 
must be derived. More important the derived eccentricity values depend on the assumed potential, 
which governs the orbit extension of a star derived from its estimated position and velocity. Another 
Achilles' heel of the method - that comparisons of eccentricities utilise a different potential for the 
calculation of orbits than the theoretical models they compare to - has in most cases been neglected. 
One may of course argue that eccentricities express the general circularity of orbits in the different 
approaches. To some part this is true, but there we come to a far more serious problem: The heating 
in models and even princip al heating mechanisms are very weakly constrained (see e.g. the discussion 
in lAumer & Binne"vll2009l) . Locally the eccentricities are quite directly related to azimuthal velocities 
in that high asymmetric drifts or respectively low rotational velocities imply large eccentricities. So 
w e see the large uncertainty i n the expected eccentricities by looking at the major changes in Fig. 2 
of ISchonrich & Binnevl d201 ll) that are induced by moderate variation in the assumed parameters of 
the disc. For example heating up the inner disc of the Galaxy a bit more, which is covered either by 
heating from a bar or respectively simply in the uncertainties intrinsic to molecular clouds as source 
of random energy, the number of high eccentricity visitors from the inner disc increases significantly. 
This would in the one-dimensional comparison be nearly indistinguishable from high eccentricity stars 
contributed e.g. by a minor merger in the outer regions of the disc. Summed up there may be physical 
entities (e.g. domination of the Solar neighbourhood by debris from some satellite) that can be ruled 
out by analysis of eccentricity distributions (though it is questionable why this should not be possible 
on full velocity space as well). However, we see that a large range of eccentricity distributions can 
be explained just by uncertainty in heating and other Galactic parameters, and moreover the use of 
eccentricity distributions in isolation hides a lot of valuable information. 



3 Conclusions 



Analysis of local spectroscopic samples reveals for the chemically young disc a significant downtrend 
of mean azimuthal velocities with metallicity, while there appears to be no notable change in veloc- 
ity dispersion. Those samples have to be used with caution as they c arry strong kinematic b iases, but 
the observations are confirmed by kinematically unbiased data from Casagr ande et alj d201 lb . While 
the qualitative picture is nicely consistent with predictions of radial migration models and has so far 
not been explained in the context of classical approaches, we stress that in particular the azimuthal 
velocity trend is not a consequence of radial migration of stars by itself. On the contrary the velocity 
trend is a natural consequence of the radial abundance gradient in the galactic disc. Radial mixing 
expands the baseline in metallicity on which the trend is observa ble, but moderates its slope. T he fact 
that the observed slope is a bit smaller than predictions from the Scho nrich & Binnevl d2009al) model 
can be reasoned by an underestimate of the amount of mixing that was caused by a too steep radial 
abundance gradient adopted in that approach. 

It is a pleasure to thank James Binney for helpful comments. 
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